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Abstract. A stability method is used to assess possible values of interspecies scattering lengths an in two- 
component Bose-Einstein condensates described within the Gross-Pitaevskii approximation. The technique, 
based on a recent stability analysis of solitonic excitations in two-component Bose-Einstein condensates, is 
applied to ninety combinations of atomic alkali pairs with given singlet and triplet intraspecies scattering 
lengths as input parameters. Results obtained for values of ai2 are in a reasonable agreement with the few 
ones available in the literature and with those obtained from a Painleve analysis of the coupled Gross- 
Pitaevskii equations. 
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1 Introduction 

Over the past decade an intensive interest is devoted to 
atomic and molecule Bose-Einstein condensation. Recent 
progress in experimental techniques makes it possible in 
the nearest future to engineer two-component Bose-Einstein 
condensates (BECs) that arc composed of two different 
alkali-metal atoms. Until now, merely theoretical stud- 
ies have been carried out for two-component alkali-metal 
atom BEC systems like Li-Rb [1], Na-Rb [1-4], and K- 
Rb [5,6]. However, the mixture Cs-Li [7,8] has been inves- 
tigated experimentally, without reaching the BEC phase. 

One of the most important quantities in engineering 
BECs is the scattering length a^ characterizing the strength 
of the interaction between the atoms of types i and j. A 
realization of two-component BECs depends on three scat- 
tering lengths, namely the intraspecies scattering lengths 
an, ^22 and the interspecies ones 012 = 021- Accurate 
knowledge of these three quantities is of great importance 
in design and implementation of two-component BECs. 
However, the problem is that the scattering lengths be- 
tween many like alkali-metal atoms are known, whereas 
those between unlike alkalis have not yet been measured 
or calculated, except for very few cases (e.g. [5-7,9-11]). 

The aim of this note is to apply a simple stability 
method to assess possible values of interspecies scatter- 
ing length by using a recent analysis [1] for existence of 
solitons within the two-component BEC. The basic idea 
of the formulation is that the solitons are stable objects 
and if they can be created in the background BEC mat- 



ter then the latter should show up also stability property. 
Various types of solitonic excitations can be created in 
two-component BEC matter, such as those of bright (B) 
and dark (D) type or of BB and DD types. Solitonic ex- 
citations in BEC matter is dealt with by many authors 
including [12-17]. In this work we shall treat only BD 
type excitations which provide finite interval for possible 
values of interspecies scattering length ai2- 

We shall determine, in tabulated form, those ranges of 
Oi2 for which various two-component BECs maintain sta- 
bility with respect to the solitonic excitations. The back- 
ground is treated in the Thomas-Fermi (TF) approxima- 
tion while the excitations are considered in the form of 
static BD solitons. Furthermore, we restrict the treat- 
ment of the dynamics to one dimension, but employ three- 
dimensional scattering lengths using an appropriate cross- 
ing area A. Note that cigar-like quasi one-dimensional 
BECs have already been produced by micro-trap devices 
developed quite recently [18-21]. As a case study we shall 
treat in detail two specific BEC systems, namely those 
composed of the alkali-metal atom pairs K-Rb and Li- 
Na. For these systems we exhibit the BD solitonic density 
profiles at the most probable values of 012 which ensure 
integrability (i.e. inverse scattering solution of the Gross- 
Pitaevksii equations), as well as the allowable range of 012 
which takes care of particle number conservation. 

The results of the present investigation indicate that 
given the intraspecies scattering length an, i = 1,2 the 
method is capable to provide possible values for the inter- 
species scattering length 012, at which the two-component 
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BEC may prove stable. Whether or not the true (physical) 
values of a\ 2 fall into these stability ranges is a question. 
If yes, then one can further calculate for finer details of 
the two-component BEC, as for the particle numbers JVj 
[i = 1, 2), the extension of components and so on. If not, 
then it would mean that other types (e.g. BB or DD) of 
solitonic excitations compatible with the physical value of 
a\2 could be realized. Moreover, one may try to adjust 
cijj's slightly by changing the strength of the magnetic 
field [22], to meet the stability requirement formulated in 
the present work. (Note, however, that the Feshbach res- 
onance method is not fully developed for interatomic col- 
lisions.) A third possibility is to choose other alkali pair 
being more favorable from point of view of values of 's 
for assembling the two-component BEC satisfying stabil- 
ity condition formulas given in the present method. 

In section [21 the theoretical background will be out- 
lined. Section 01 contains the estimated ranges of inter- 
species scattering length a\ 2 , in tabulated form, for which 
the solitonic stability requirement holds and a compari- 
son with earlier findings and those resulted from Painleve 
analysis will be made. SectionQ]is devoted to a short sum- 
mary of the method. 



2 Theoretical background 

Following [1] we write the coupled Gross-Pitaevskii (GP) 
equations [23,24] describing a quasi one-dimensional two- 
component BEC in the form (i = 1,2): 



2m, 
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has been applied earlier by Dum et al [12].) Inserting the 
above ansatz into the GP equations Q and assuming that 
the excitation mechanism restricts the change of density 
only to a small interval around x — where the TF solu- 
tions can be approximated by <Pi(0), the following coupled 
equations can be obtained for the perturbing functions 
(t = l,2) 



where denotes the individual mass of the ith atomic 
species, Qij = 2ith 2 'aij / 'A^bij with aij being the 3D scat- 
tering length, A represents a general transverse crossing 
area of the cigar-shape BEC, fj,ij = miirij / (m, + mj) 
is the reduced mass, and Vi, (i = 1,2) denotes the ex- 
ternal trapping potentials. In the case of real trap po- 
tentials the normalization of the wave functions reads as 
Ni — \ ipifdx, (i = 1, 2) with JV, denoting the number 
of atoms in the ith component. We do not allow for the 
species to transform into each other, so that the numbers 
of particles are conserved quantities for both alkalies. 

The stationary solutions of the coupled GP equation 
Q can be written in the form^i (a;, t) = <&i(x) exp(— iEit/K) 
where Ei represents the one-particle energy of the i— th 
component particle [i = 1,2). Using this formula and ne- 
glecting the kinetic terms one is able to derive the approx- 
imate density profiles and a semi-infinite range for a\ 2 in 
which the Thomas-Fermi ground state exists. In the case 
of a harmonic trap potential the calculation results in the 
usual inverse-parabolic distribution. 

We then look for the excited static solutions in the 
form ipi(x,t) = 4>i(x)4>i(x) exp(—iEit / H) with <fii(x) being 
an excess or defect of the i— th component of the back- 
ground density, triggered by an excitation mechanism. 
(Note that such trial function for a solitonic excitation 
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1,2) and the very small potential 
terms have been neglected. 

The coupled equations above can be solved for the per- 
turbing functions 4>i(x) by using the static bright-dark 
(BD) coupled solitonic ansatz 



i>Bi{x) = <?isech(fca;), 4>bi(x 
i>D2 (x) = Q2 tanh (kx) , 0n 2 (x 



±oo) = 0, (3a) 
±oo) = ±q 2 , (3b) 



with the specified boundary conditions, generic complex 
amplitudes qi, and common real size parameter k. 

The requirement that the moduli of the two amplitudes 
qi and q 2 be positive real numbers gives the stability con- 
ditions (A^ = dij(\ + mi /raj)) 
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for the existence of BD solitonic excitation within the two- 
component BEC. The finite ranges of a\ 2 provided by the 
solitonic stability requirement (|4al -b) are naturally nar- 
rower than the TF domain of a\ 2 prescribed by the posi- 
tivity of the density | <Pi (x) | 2 . 

Let us note here that instead of the BD solitonic ansatz 
(|3al -b) one might as well apply BB or DD prescriptions 
which would still satisfy equations J3J). These ansatzes re- 
sult in similar stability conditions like the forms (|4al -b) but 
yield different stability domain of a\ 2 . This is one reason 
for the usage of the terminology assessment by which we 
mean that the physical value of a\ 2 may be well outside 
the ranges given later in this paper. Nevertheless, for the 
specific conditions corresponding to the BD stability we 
shall find compatibility between our findings and others, 
including the P-test a\ 2 values which are obtained from 
the singularity structure analysis of the coupled GP equa- 
tions. 

From the many possibilities we shall select two exam- 
ples, with alkali pairs of Li-Na and K-Rb, representing 
realistic possibilities for future experimental study. 



2.1 Static BD solitonic excitations of Li-Na and K-Rb 
condensates 

In order to illustrate the solitonic excitations we show in 
figure the density profiles ^(a;)! 2 of the two BEC sys- 
tems composed of lithium-sodium and potassium-rubidium 
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atoms. All of the components are chosen to be in the 
triplet scattering state. We have kept N2 fixed by choos- 
ing the dark soliton particle number to be N2 — 10 4 for 
both the 7 Li- 23 Na triplet system and the 41 K- 87 Rb two- 
component condensate. The bright component has parti- 
cle number Ni rs 5 x N2 in both cases, corresponding to 
the chosen values of an = 2.9 nm (Li-Na) and 012 = 4.85 
nm (K-Rb). The excitation of the BEC is manifested by 
the appearance of the B and D solitons on the TF back- 
ground density. 



7 Li- 23 Na system, and from ayi — 4.1; 7.5 nm to ayi = 4.5; 
7.4 nm for the 41 K- 87 Rb configuration. 

In conclusion we can say that if the actual value of a\2 
determined by some experimental or theoretical method 
turns out to fall within the range determined by inequal- 
ities <|4al -b) then one may fine tune the particle number 
ratio N2/N1 and the size parameter k. In this way we 
get complete information about both the background and 
solitonic stability. 



2.2 (7 12 values resulting from Painleve analysis 
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Fig. 1. Absolute squares of wave functions, l^i] 2 an d \ip2 1 2 in 
unit [fim -1 ], of the BD solution of the coupled Gross-Pitaevskii 
equation as a function of the position x for different BEC 
systems, (a) triplet 7 Li- 23 Na system: N 2 = 10 4 , Ni » 4.7 x 10 4 , 
an = -1.4 nm ( 7 Li), a 2 2 = 4.0 nm ( 23 Na) ai 2 = 2.90 nm, trap 
frequencies are: uji = 2n x 530 Hz, L02 = 2n x 292 Hz. (b) triplet 
41 K- 87 Rb system: N 2 = 10 4 , N x as 4.9 x 10 4 , a 1± = 3.4 nm 
( 41 K), 022 = 5.5 nm ( 87 Rb), a±2 = 4.85 nm, trap frequencies 
are: ui = 2n x 220 Hz, uj 2 = 2n x 150 Hz. Figures (a2) and (b2) 
show the dependence of the logarithmic ratio ln(A^/iVi) on 
the interspecies scattering length a\2 for these two-component 
BECs. The transverse crossing area A — r 2 iv with radius r = 
0.3 nm has been employed throughout. 



If the system parameters mi, 7772, an and 022 are 
given, then one can determine the broadest range of the 
interspecies scattering length by using the inequalities (|4al 
b). This procedure does not involve the particle numbers 
A^i and N2 and gives an assessment of ai2 range for which 
the existence of BD or DB solitons in the two-component 
BEC might be expected. Particle number conservation fur- 
ther restricts the interval of 012 and leads to a system cou- 
pled equations for N2/N1 and the size parameter k (see 
eqs. (11) of Ref. [1]). Solving these equations, the final re- 
sults are shown in figures^a2) and^b2). We observe here 
a rather smooth behavior of N2/N1, as well as a shrink- 
ing of the allowable intervals from a\2 = —0.6; 6.1 nm (as 
determined by inequalities (@a}b)) to a\2 — 1.9; 4.9 nm 
(obtained by use of particle number conservation) for the 



In Ref. [25] it has been shown that the coupled GP equa- 
tions pass the Painleve (P) test provided the following 
condition among the system parameters (masses, interac- 
tion strengths) is satisfied 



(mifln + m 2 a 2 2) . 

012 = ; ; — ± 

M{m 1 +m 2 ) 



(5) 



with M 



(mian + m 2 a 2 2) 2 + 4M(M - 2)mim 2 an022 
M(mi -I- m 2 ) 

(2m + l) 2 + 7 /16, where m can be treated 



as an arbitrary positive integer number (m = 0, 1, 2, . . .). 
This expression depends only on the ratios 7711/777-2, an/ai2i 
and 012/022 involving the characteristic parameters of the 
BEC, and the external potential families can be classi- 
fied by the different values of m. Experimentally preferred 
spatially harmonic (~ x 2 ) potential- families fall into the 
category 777 = 2 (M = 2). The condition given by equa- 
tion (JSJ can be written in the following interesting form 
M = A 2 2 7,02(012) — ^4ii/si(oi2)- This relation connects 
two independent analysis, namely the solitonic stability 
conditions (0aj-b) and the P test encoded into equation JSJ) . 
The left hand side is responsible for the fulfillment of the 
P test which constitutes one of the ingredients of integra- 
bility conditions [26] meaning among many others that 
solitonic solution of the GP equations can be found by 
the method of inverse scattering transform. Right hand 
side contains the functions fsi and jui which should be 
positive in order that the coupled GP equations may have 
static solitonic solutions of BD or DB types. 



© 3 Results and discussion 

3.1 <M2 values allowing stability and BD solitonic 
excitations of BECs 

In this subsection the BD stability ranges of 012 will be 
presented in tabulated form, using both triplet and singlet 
intraspecies data as input, for ninety pairs of alkali atoms. 
Note that the constraint arising from TF stability is not 
included in the data given. 

TablesHandElgive, separated from each other by semi- 
colons, the limiting point values of the intervals of the 
interspecies scattering length 012 satisfying the stability 
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conditions encoded into inequalities (|4al -b). These condi- 
tions express the possibility of developing solitonic excita- 
tions of static BD or DB type from the ground state of the 
two-component BECs, irrespective of the particle num- 
bers Ni and N2 of the mixture. The data have been cal- 
culated, respectively, with triplet and singlet intraspecies 
scattering lengths an (i = 1,2) listed in the columns 
headed by the labels at and a s . The exception is the 
133 Cs atom where only the average scattering length (an — 
+ ^/cr e xp/8Tr) has been available (therefore the parenthe- 
sis). 

In general we observe in tables ^ and [21 that one of the 
limiting points of the intervals of the corresponding DB 
and BD cases coincides; a fact easily derivable by using the 
relation = /_Bi(ai2) = /rata^)- The diagonal elements 
of tables ^ and [21 represent one-component BECs which 
is out of scope of the present interest so that the data 
are missing there. Furthermore, the rows corresponding to 
7 Li, 39 K, and 85 Rb of table Q do not contain any results 
as well, if their partner owns a positive triplet intraspecies 
scattering length. The lack of information means here that 
we do not find any region of triplet interspecies scattering 
length ai2 which gives rise to a stable DB type excitation. 
On the other hand all items of the column headed by 7 Li, 
39 K, and 85 Rb are filled by data. This means that there 
are definite ranges of triplet a\2 which support BECs with 
excitations of BD type where the 7 Li, 39 K, and 85 Rb atoms 
play the role of the B component (owing to their attractive 
interaction character). The above argumentation can be 
confirmed by using stability conditions (^ajb) with an < 
and 022 > 0. 

In the case when both an and 022 are positive we ob- 
serve the existence of BECs with both DB or BD forma- 
tions within a finite interval of positive interspecies scat- 
tering length ai2- This can be understood physically in 
the following way. In the context of the mean field theory, 
an effective potential is created by the i-th atoms acting 
repulsively among each other (positive scattering lengths) 
and thus forming the cavity of the D component. Inside 
this cavity there is a possibility for the j-th atoms to form 
a bunch representing the B component despite their in- 
herent repulsive character. 

The values of 012 given in tables ^ and [21 can be used 
for orientation and in design of two-component BECs. 
One starts from known intraspecies scattering lengths an 
(i = 1, 2) which are represented by the data denoted by at 
or a s in the tables. Then one performs the range determi- 
nation by using inequalities (|4"al b). Now one is provided 
with guessed values of interspecies scattering length ai2 
which can maintain stable two-component BEC configu- 
ration. If, from some other sources, we are assured that the 
actual physical values of ai2 falls outside the range calcu- 
lated then one may try either to move the a^- values into 
the proper direction by utilizing the Feshbach-resonance 
procedure [22] or to carry out BB or DD stability analysis 
to learn what type (if any) of static solitonic excitations 
of the components are possible (suggesting, at the same 
time, stability of the ground state configurations as well) . 



However, if all a^ 's = 1, 2) satisfy inequalities H4a^ 
b), then one may start to carry out the production of 
two-component BEC using information obtained by fine 
tuning the component numbers Ni as described in Rcf. [1]. 

3.2 ai2 from earlier studies 

There are very few interspecies scattering length ai2 data 
in the literature. The cases studied so far include the pairs 
composed of different isotopes of rubidium atoms [10], of 
those of potassium atoms [9] , of potassium and rubidium 
atoms [5,6], of cesium and rubidium atoms [11], and of 
lithium and cesium atoms [7]. In the latter case the sign 
of ai2 remains undetermined since we can derive, by using 
the relation a e \ = Aira\ 2 , only the absolute value of the 
interspecies scattering length from the measured unpolar- 
ized elastic cross-section data a e p( 7 Li — 133 Cs) « 5 • 10~ 12 
cm 2 to be ai 2 ( 7 Li — 133 Cs) = ±6.3 nm. This value of 
ai2 is to be compared to our ranges of ai2 character- 
ized by the limiting points —0.1; 4.6 and 0.9; 4.6 nm ob- 
tained, respectively, with triplet and singlet intraspecies 
parameters an = —1.4 and 1.7 nm for the lithium atoms, 
and 022 = 2.4 nm for the cesium atom (where the latter 
value also represents unpolarized data [27] with known 
positive sign). From the comparison we predict the plus 
sign for the unpolarized interspecies scattering length, i.e., 
ai2( 7 Li — 133 Cs) = +6.3 nm, provided the stability can be 
attained. 

For the potassium-rubidium system we compare our 
ranges ai 2 ( 39 K - 87 Rb) = -0.6; 7.6 nm and ai 2 ( 41 K - 87 
Rb) = 4.1; 7.5 nm obtained using triplet intraspecies data 
(i.e. with an( 39 K) = -0.9 nm, an( 41 K) = 3.4 nm, and 
fl22( 87 Rb) = 5.5 nm) to the triplet interspecies values 
ai 2 ( 39 K - 87 Rb) = 1.6 nm [5] and a i2 ( 41 K - 87 Rb) = 
8.1J3' 5 nm given in Ref. [6]. We see here that our findings 
are comparable with the calculated values of ai2- 

For the cesium- rubidium system a recent work [11] lists 
various ai2 values calculated using six different ab-initio 
potential for both triplet and singlet collisions. The cal- 
culated data lying between 012 = 2; 19 nm for the singlet 
and between ai2 = —8.5; 3.2 nm for the triplet 133 Cs- 85 Rb 
system compare well to ours given by the BD ranges of 
a 12 ( 133 Cs - 85 Rb) = 16.9; 97 nm and ai 2 ( 85 Rb - 133 Cs) = 
3.0; 16.9 nm for the singlet and ai 2 ( 85 Rb - 133 Cs) = 
— 14.8; 2.9 nm for the triplet arrangement, respectively. 
In the case of the 87 Rb- 133 Cs DB system the singlet scat- 
tering lengths are comparable; for example, the calculated 
singlet value of ai2 = 3.0 nm of the ab-initio result VI [11] 
compares well to our range of ai2( 87 Rb — 133 Cs) = 2.9; 3.3 
nm. 

Finally we mention the singlet interspecies scatter- 
ing length for different isotopes of the potassium atoms, 
namely the value ai2( 39 K— 41 Kb) = 5.9 nm [9] to be com- 
pared to our range ai2( 39 K — 41 K) = 5.7; 7.0 nm of the 
DB case. Again, a clear compatibility is observed. 

We may thus conclude that the solitonic excitation 
method carried out in section [21 is capable to assess ap- 
proximate intervals of interspecies scattering length ai2 
which may encompass the actual values themselves. 
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3.3 ai2 values derived from Painleve test 

In table the values of an calculated from equation iJSJ 
are shown for the harmonic oscillator potential obtained 
with singlet (upper diagonal) and triplet (lower diago- 
nal) intraspecies scattering lengths an, i = 1,2. Although 
ai2 = is always a solution of equation at M = 2, 
this value is not shown in table |3| because it represents 
the uncoupled case. 

Since the P test of a nonlinear partial differential equa- 
tion primarily serves to explore its singularity structure, it 
is only implicitly connected with soliton formation. There- 
fore we do not in general expect that all values of table 
match with those encompassed by the a\2 intervals of ta- 
ble and |21 Yet, the values of items in tables and 
(lower diagonal) and t ablest and (upper diagonal) com- 
pare well. From this comparison we can learn that only 
those DB excitations are possible in which the D compo- 
nent is represented by the heavier elements listed in the 
first line of table [3J Exceptions are the 7 Li atom which 
always plays the role of the B component, and the 85 Rb 
atom which always constitutes the role of the D compo- 
nent (when singlet scattering is considered) . This result is 
physically reasonable, since in a mixture of atoms with dif- 
ferent masses and comparable (positive) interspecies scat- 
tering lengths, the atoms of the slower (more massive) 
component is likely to form a dip of D soliton, into which 
the atoms of the swifter (lighter) component can be cap- 
tured to create a B soliton. The exceptions occur if the 
repulsion/ attraction between like atoms is strong enough 
as in the case of 85 Rb/ 7 Li to cause a situation with per- 
manent D/B components. 

4 Summary 

The assessment of the interspecies interaction parameter 
<ii2 at which two-component atomic BEC exhibits stable 
configuration is an important task, and can be carried out 
with the present method in the following way. 

In principle, all the ground state parameters includ- 
ing a\2 can be extracted from the TF equilibrium den- 
sity profiles. This method is rather ambiguous and, there- 
fore, restricting the number of uncertain parameters or 
their ranges may prove useful in design of two-component 
BECs. Moreover, in this fitting procedure the particle num- 
bers are not coupled, they can be chosen independently 
of each other. The coupling is represented mostly by the 
inter- and intraspecies scattering length parameters 
(i, j = 1, 2) being strongly correlated with each other. 

If one is interested in producing excitations of BEC in a 
form of static solitons then one first determines a region of 
interspecies scattering length ai2 which may support the 
soliton formation. This task can be performed by simulta- 
neous solution of the inequalities (TJaJ-b) for 0,12, provided 
the interaction parameters an and (122 are given. This step 
is again free of the particle number parameters iVi and N2 
and gives the broadest region of a±2 for which DB or BD 
static solitons ever can be created inside the stable two- 
component BEC. For the two examples treated explicitly, 



table n contains such intervals between a\2 = —0.6 : 6.1 
nm for the 7 Li- 23 Na triplet system, and 0,12 = 4.1; 7.5 nm 
in the case of the 41 K- 87 Rb triplet scattering. 

If the actual value of intraspecies scattering length 012 
which must be known from other source, is within the in- 
terval determined above, then one may perform a fine tun- 
ing for the particle number ratio N2 /N\ at which the soli- 
tonic excitation may be implemented. This fine tuning is 
accomplished by solving the coupled nonlinear equations 
(11) of Reference [1] for the ratio N2/N1 and the size pa- 
rameter k of the solitons. For the two explicit examples 
chosen, figures ^ contain the results with N2 kept fixed. 

Bright-dark soliton combinations in one- and two-comp- 
onent condensates have already been considered in a vari- 
ety of papers, including [12, 14, 15, 17]. Experimental cre- 
ation of both bright [16] and dark [13] solitons have been 
observed in case of one-component condensates. We hope 
that fine tuning the system parameters (scattering lengths, 
particle numbers, trapping frequencies) as outlined above 
may also help in realizing coupled solitonic excitations in 
two-component Bose-Einstein condensates. 

This work has been supported by OTKA under Grant Nos. 
T038191, T047035 and M36482. 
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Table 2. Intervals of singlet interspecies scattering lengths 012 assessed by static solitonic excitations of two component BECs composed of alkali pairs indicated. 
The column as gives the singlet scattering length an of atoms taken from the references: 7 Li [29,32,33,41], 23 Na [34,35], 39 K, 41 K [9], 83 Rb, 85 Rb, 87 Rb [10], 133 Cs, 
135 Cs [27,39]. (For X H the triplet value has been used [28]). All values are in unit of nm. 
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Table 1. Intervals of triplet interspecies scattering lengths ai2 assessed by static solitonic excitations of two component BECs composed of alkalis indicated. The 
column a t gives the triplet intraspecies scattering length a u taken from the references: X H [28] , 7 Li [29-33], 23 Na [34-36], 39 K, 41 K [9,37], 83 Rb, 85 Rb, 87 Rb [10,38], 
133 Cs, 135 Cs [27,39,40]. All values are in unit of nm. 
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Table 3. Values of interspecies scattering length a\2 solving the P test condition, equation JjjJ, at classification numbers m = 2 preferring harmonic oscillator trapping 
potentials. Upper diagonal contains results obtained with singlet intraspecies scattering lengths as listed in the second column, lower diagonal exhibits those with 
triplet da's given in the third column as a^. All values are in unit of nm. 
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